Orbital evolution of a particle around a black hole: II. Comparison of contributions of 

spin-orbit coupling and the self force 
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We consider the evolution of the orbit of a spinning compact object in a quasi-circular, planar orbit 
around a Schwarzschild black hole in the extreme mass ratio limit. We compare the contributions 
to the orbital evolution of both spin-orbit coupling and the local self force. Making assumptions 
on the behavior of the forces, we suggest that the decay of the orbit is dominated by radiation 
reaction, and that the conservative effect is typically dominated by the spin force. We propose that 
a reasonable approximation for the gravitational waveform can be obtained by ignoring the local 
self force, for adjusted values of the parameters of the system. We argue that this approximation 
will only introduce small errors in the astronomical determination of these parameters. 
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I. INTRODUCTION AND SUMMARY 

One of the most promising sources of gravitational ra- 
diation that can be detected by the Laser Interferometer 
Space Antenna (LISA) Q is the capture of a compact ob- 
ject by a supermassive black hole. In the extreme mass- 
ratio limit this problem can be solved using perturbation 
theory, in which the compact object (e.g., a stellar mass 
black hole or a neutron star) is treated as a test particle 
moving on the fixed background of the central, super- 
massive black hole. 

A test particle, in the absence of self interaction and in- 
ternal structure, moves along a geodesic of the spacetime 
background. When the particle's mass is much smaller 
than the typical length scale of the curvature of the back- 
ground spacetime, it is useful to use perturbation theory. 
The particle's mass (or charge, if any) contributes to the 
spacetime curvature, such that the particle moves along a 
geodesic of a perturbed spacetime. (These perturbation 
fields can be interpreted as those of free gravitational 
waves, which are produced by the particle at retarded 
times. These perturbations are smooth on the particle's 
worldline 2].) Alternatively, the particle can be con- 
strued as moving along an accelerated, nongeodesic tra- 
jectory of the unperturbed background spacetime: Let 
the mass of the companion compact object be \i, and the 
mass of the central black hole be M, such that /i <C M. 
The particle's acceleration at order \x is driven by forces 
at order /x 2 imparted on the particle. To obtain the cor- 
rected orbit of the particle [to 0(/i)] one needs to include 
all the forces, at 0(/i 2 ), which act on the particle. 

Two sources can contribute forces at 0(/i 2 ): i) the 
particle's self force Q, and ii) the particle's spin angular 
momentum: When the companion of a central black hole 
is a black hole of mass fi, the latter carries spin angu- 
lar momentum j = s/i 2 , where — 1 < s < 1. (Neutron 
stars also have maximal spin angular momentum which 
is quadratic in their mass J4|.) The force imparted on the 
particle is at order /i 2 , which endows the particle accel- 
eration at order [i. Numerous authors have considered 
the orbital evolution under either self force effects [j| or 



spin effects 6]. However, in order to obtain the parti- 
cle's orbit to order \x, one needs to include both the self 
interaction and the spin effects. 

The relative importance of the two contributing effects 
can be evaluated using an ordcr-of magnitude estimate 
for r 3> M. The component of the self-force responsible 
for the decay of the orbit (dissipation, "radiation reac- 
tion") scales like (fi/M) 2 (M/r) 5 7}. The contribution 
of the particle's spin angular momentum to the decay 
of the orbit can be found from the Papapetrou equa- 
tion. Specifically, an order-of magnitude estimate for the 
spin force can be found from the fact that it is linear in 
the particle's spin, linear in the spacetime curvature, and 
quadratic in the particle's four-velocity. (One needs to be 
careful here, because the various components of the four- 
velocity scale differently in M/r. Specifically, u f ~ 1, 
whereas u v ~ (M/r 3 ) 1 / 2 — see below.) Taking the typi- 
cal curvature to scale as M/r 3 , and the typical velocity of 
a particle in a quasi-circular orbit to be given by Kepler's 
law, one expects the contribution of the spin force to dis- 
sipation to scale as (^t/M) 3 (M/r) 13 / 2 . Consequently, we 
expect the decay of the orbit to be dominated by the ra- 
diation reaction, with only small corrections (which we 
neglect) due to the particle's spin angular momentum. 

The situation is much different when one considers the 
conservative correction to the orbit at 0(/x). The com- 
ponent of the self force contributing to the conservative 
effect scales like (/x /M) 2 (M/r) 6 , whereas the component 
of the spin force contributing to the conservative effect 
scales like (/z/M) 2 (M/r) 7 / 2 . (See below the justification 
for these scaling laws.) Both forces are at 0(/i 2 ), but be- 
cause of the slower drop off with M/r, naively one might 
expect the spin force under these conditions to be much 
more important than the self force. In a strict sense, how- 
ever, a direct comparison of the spin force and the self 
force tells us only little about the relative importance of 
the two forces. The reason is that the self force is gauge- 
dependent. In particular, one can always choose a gauge 
in which the self force vanishes. Consequently, the scal- 
ing law used above is not unambiguous, at least as long 
as we do not specify the gauge in which it is written. 
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A useful way to compare the spin force and the self 
force then is to compare between gauge-independent 
quantities. The generated waveform or the number of 
cycles that the system spends in a logarithmic interval of 
frequency are two such quantities. In practice we evolve 
a quasi-circular, equatorial orbit around a Schwarzschild 
black hole using both forces, or either one, and compare 
these gauge-independent quantities. [An initially equato- 
rial orbit will remain nearly equatorial because the total 
angular momentum is conserved, and the spin angular 
momentum is much smaller than the orbital angular mo- 
mentum. Specifically, the spin angular momentum is at 
0(/z 2 ), such that its rate of change is at 0(p?). Conser- 
vation of total angular momentum imply that changes in 
the orbital angular momentum are also at 0(/i 3 ), whereas 
the orbital angular momentum itself is at 0([a). Conse- 
quently, changes in the direction of the spin vector imply 
only very small changes in the direction of the orbital 
angular momentum, or very small changes in the orbital 
plane §.] Indeed, we find that our naive expectations 
are realized: the conservative effect is controlled by the 
spin force, which overwhelms the self force. More general 
orbits around a Kerr black hole will include also dissipa- 
tion due to the spin force and spin-spin coupling. Here, 
we ignore such effects, and focus on the spin-orbit cou- 
pling. Our motivation is that for such a simplified orbit 
we can already demonstrate our main point: spin forces 
for astrophysical systems may be much more important 
than the self forces, as far as conservative effects on the 
orbital evolution are concerned. 

In particular, we find that neglecting the local, con- 
servative self force in the construction of templates may 
introduce only a small error in the astronomical determi- 
nation of the parameters of the system. Specifically, we 
can match a template made without the local self force 
almost exactly to a waveform which does include the lo- 
cal, conservative self force, but with a slightly different 
value for the spin of the companion. 

In a more general system, the component of the spin 
force which induces the conservative effect also con- 
tributes to dissipation. Consider a Kerr black hole in 
Boyer-Lindquist coordinates. The rates of change of the 
particle's energy, E, and angular momentum (in the di- 
rection of the spin axis), L z , can be found from the fluxes 
of energy and angular momentum to infinity and through 
the central black hole's event horizon. The rate of change 
of Carter's constant, Q, however, cannot be found us- 
ing balance arguments and global conservations laws, be- 
cause Q is non-additive. To find Q one needs to find first 
the local forces which act on the particle. Specifically, 

v r 

Q = G E (COM,g)E + G Lz (COM,g)L z - 2£— f r , 

W 

where Ge, Gl, are certain (known) functions of the con- 
stants of motion (COM) in the absence of dissipation, 
and the metric g and S = r 2 + a 2 cos 2 9, a being the 
spin parameter of the black hole. To have a full descrip- 
tion of dissipation then, one needs the total radial force, 



f r , which acts on the particle. Similarly to our argument 
above, we expect the contribution to Q due to the self 
force to be small compared with the contribution of the 
spin force. 

As generically the companion is expected to spin (and 
even spin fast), we propose that a reasonable approxi- 
mation for the orbital evolution of a spinning particle 
can be found without finding first the local self force. 
Specifically, we propose that a practical way to obtain 
the orbital evolution (and the waveforms) to high accu- 
racy could be the following: First, find E, L z using the 
fluxes to infinity and down the event horizon of the cen- 
tral black hole. Next, find Q using E, L z and the spin 
contributions to f r . Undoubtedly, that will introduce an 
error in the determination of Q. However, because of 
the smallness of the radial component of the self force 
with respect to the radial component of the spin force, 
we believe that this may be a useful approximation for 
the total Q. To find the orbital evolution one needs also 
to include conservative effects [13 . Again, we propose to 
neglect the contribution of the self force to the conserva- 
tive effects, and approximate the full conservative effect 
by the conservative effect due to the spin force. 

In the simple problem we consider here the motion 
is quasi-circular and equatorial around a Schwarzschild 
black hole, such that the spin angular momentum of the 
companion is aligned (or anti-aligned) with the orbital 
angular momentum. It is not hard to let the central 
black hole be spinning — with the spin axis pointing 
along the same direction as the companion's spin angular 
momentum and the orbital angular momentum — and 
include also spin-spin coupling effects. We are hoping 
to return to that problem in the future. We note, that 
when the spins are not aligned as in our assumptions, the 
precession of the particle's spin will induce a changing 
quadrupole moment, such that the spin force will cause 
a strong dissipative effect. The implications of such a 
spin dissipative effect are beyond the scope of the present 
paper. 

Our goal in this paper is more to point out that the ac- 
curate determination of the local self force may perhaps 
be less crucial for determination of the orbital evolution 
of extreme mass-ratio binaries than previously thought, 
than to give a definite prediction. First, some compan- 
ions may have only little spin angular momentum, in 
which case one would no longer be justified in neglect- 
ing the local self force. Second, and most importantly, 
our analysis in what follows is incomplete, in the sense 
that some of the necessary pieces of information for car- 
rying out the comparison of the spin and self force effects 
are as yet unknown. We consequently make a number of 
assumptions, that allow us to obtain what we believe to 
be at least a reasonable order-of-magnitude estimate for 
the effect of interest. We emphasize that our analysis is 
based on a number of assumptions, which are as yet un- 
proven. We believe that our main conclusion — namely, 
that spin-orbit coupling may overwhelm the conservative 
self force effect — is insensitive to the accuracy of these 
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assumptions and is robust. At three places our assump- 
tions are speculative (see below for details): First, we as- 
sume that the terms linear in second-order forces in the 
expression for the second-order radial velocity are small 
compared with terms quadratic in first-order forces. Sec- 
ond, we assume that finite-mass effects in the luminosity 
in gravitational waves, which are small in the weak field 
regime, are small also in the strong field regime. Third, 
we assume that the radial component of the gravitational 
self force is proportional to its scalar field counterpart, 
with a specific proportionality constant. Of these three 
assumptions, the last one is the easiest to test. When this 
is done, the actual gravitational self force can be used to 
replace the assumed expression l|27|l. While these as- 
sumptions appear to be quite natural to make, they are 
by no means guaranteed to be justified under all circum- 
stances. However, even if any of them turns out to be in- 
correct, the implications on our main point in this paper 
are not expected to be strong. We therefore believe that 
our results are at least reasonable order-of-magnitude es- 
timates, and that our main point is relevant for realistic 
extreme mass-ratio binaries. 

The organization of this paper is as follows. In Sec- 
tion ^ we describe the equations of motion and derive 
the equations for a quasi-circular, equatorial orbit with 
aligned spins, using the Papapetrou equations and lin- 
earizing in the spin covector. Appendix 1X1 includes more 
details on the definitions of the spin force and the local 
self force. In section ITTT1 we derive the perturbative solu- 
tion for the equations of motion following the method of 
Ref. (paper I), but with slightly different definitions 
and notation, which will make the generalization to a 
Kerr black hole simple. Then, we discuss the domain 
of validity of the solution. In Section IIVI we describe 
the computation of the evolved orbit and the number of 
cycles that the system spends in a logarithmic interval 
of frequency. Next, in Section El we discuss the model- 
ing of the self force. First, we discuss the fitting of the 
numerically-derived luminosity in gravitational waves to 
a smooth function using a match to two different asymp- 
totic expansions, and discuss the associated errors, and 
then we discuss our conjecture for the as-yet undeter- 
mined (radial component of the) self force. Finally, in 
Section IVII we compute the orbit and the waveforms, 
and discuss the relative importance of the spin force and 
the self force. 



II. EQUATIONS OF MOTION 
A. Equations of motion to order /i 2 



larly, the deviation of the orbit from the geodesic because 
of the self force is also at order fi, such that its effect on 
the spin force is again negligible. It is clear then, that the 
total force, at order /i 2 , which acts on the particle and 
pushes it off the geodesic is just the sum of the self force 
and the spin force. Specifically, the equations of motion 
for a particle of mass /i, whose center of mass travels 
along the worldline with four- velocity it M = dz^ / dr 
are given by 



/'- 



Du a 
~drT 



J = /SF /sDi 



/ spin ' 



(i) 



The expressions we use for /c? F and /^; n are given in 
Appendix [X] 



B. Specializing to equatorial motion with an 
aligned spin 

1. General equatorial motion 

The line element, in the usual Schwarzschild coordi- 
nates, is given by 



ds 2 



-F(r) dt 2 



dr 2 

no 



] dn 2 



(2) 



where F(r) = 1 - 2M/r, and dil 2 = d6 2 + sin 2 9 dip 2 . 
(Notice, that the form of the Papapetrou equations of 
motion depends on the choice of signature.) 

In general, the spin covector precesses along the mo- 
tion, and the orbit is not planar. In particular, with 
arbitrary alignment of the spin covector, an otherwise 
equatorial motion will no longer be equatorial. To facili- 
tate the analysis, let us specialize to a particular solution 
of the Papapetrou equations, in which the entire motion 
is on the equatorial plane. Under the requirement that 
the entire motion is equatorial, Eq. (| A9|) for the paral- 
lel transport of the spin covector becomes very simple. 
Specifically, Eq. I|A9|I becomes 



dS a 
~dT 



1 



u r S a St 



(3) 



which implies that S e = sfi 2 /r, and 5* = S r = S* = 0, s 
being a constant. Notice that the magnitude of the spin 
vector then is S = sfj, 2 , such that s satisfies — 1 < s < 1 



(4) 



for black holes. Equation l|A8(l then becomes 



~dr 



M 



F(r 



■u* u r S e 



The total force which acts on the particle at order /x 2 
is just the sum of the well known spin force and the self 
force. The deviation of the orbit from the geodesic be- 
cause of the spin force is at order [i. Its effect on the self 
force is therefore at order /x 3 , and hence negligible, as we 
are interested here in the forces only at order fj 2 . Simi- 



Du r M t e 
— - — = 3 — b (r)u u v S 
dr r 



/'- 



dr 



(5) 



(6) 
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FIG. 1: A black hole with mass fi, spin angular momentum 
j, and tangential velocity v in circular, equatorial motion 
around a black hole of mass M. The orbital angular momen- 
tum is either aligned or anti-aligned with the spin angular 
momentum j . In this figures the angular momenta are shown 
aligned. 



Di 



dr 



(7) 



Notice that the last equation indeed satisfies the require- 
ment for equatorial motion. The equations of motion 
|@J can be interpreted as an external force, which is 
given by their RHS, acting on the particle and pushing 
it off the geodesic. This spin force joins the self force, as 
they are both at order /i 2 . 



2. Quasi-circular equatorial orbits 

The equations of motion simplify even further for the 
case of quasi-circular orbits. For such orbits the radial 
motion is perturbatively small, i.e., the radial velocity is 
at order /j,. That is the case when an initially circular 
orbit damps under radiation-reaction effects. Keeping 
terms on the RHS of Eqs. Q-Q to order [i 2 , we find 
that the only remaining component of the spin force is 
the radial component. The spin force we use then is given 
by 



ret 
J spii 



M 



2M 



u t u v S e 6? + 0(m 3 ) . (8) 



III. PERTURB ATIVE SOLUTION 

A. Derivation of the solution 

We next solve the equations of motion (JIJ pertur- 
batively, for the case of quasi-circular, equatorial orbit 
We use the normalization condition for u a , namely 
u a u a = —1, to eliminate it* from the equations of mo- 
tion. We next use the t component of the equations 
of motion (EOM) to eliminate We can simplify the 
EOM to first-order (nonlinear) ODEs by taking r — V(r), 
x = Vx'(r), and x denotes any quantity. We find the 
EOM to be 



VV' - 



V 2 

3M — + (M - Lu 2 r 3 ) 



A 

7i 



1 



A 2 V 
jiu \r 2 A 



VJ + 2. { r-ml= {r - 2M) {^ 
A fiu 1 Ar 



(9) 



(10) 



where u> is the angular velocity, A 
,t 2 — 1 _ o/w/, _ „T/2 



r 2 - 2Mr, and 

l/u tA = 1 - 2M/r - rV 2 /A - r 2 u 2 . We next expand 
Eqs. ltTU|) in powers of e = fi/M: u = W( ) + Wm + 

w (2), V = V(i) + V(2), and a l = a-p + af\ Here, X( 3 ) 
denotes the term in x which is at 0{e J ), and a.; being 
the self acceleration. We then expand the self force as 

fF = fi 1] + fP, where fi j) = We then expand 

Eqs. (f^Tfl — (fTOf) . and solve perturbatively order by order. 

The zeroth order term of Eq. © recovers Kepler's law. 
Specifically, it yields 



J (o) 



Mi 

7"' 2 



(11) 



The first order corrections to Kepler's law are obtained 
from the terms at 0(e) of Eq. Specifically, we find 
that 



(i) _ 



rfr 



(1) 



(r - 2M) 



2/y 



r 3 ^(o) 



(12) 



We next define 

vj 2 := r 2 Au[ 0) + 2r 2 (r - 3M)uj (0) . (13) 

The first order term of the radial velocity is obtained 
from the terms at O(e) of Eq. lfTU|) : 



We thus find that to order fj, 2 , the only contribution to 
the temporal component of the total force comes from the 
self force: dissipation is controlled by radiation reaction. 
The conservative effects, however, have contributions (via 
the radial force) from both the self force and the spin 
force [as they are both at 0(fi )]. In fact, we find that 
the spin contribution to the radial force is typically much 
larger than the self force contribution. We shall study 
their relative importance in detail below. 



V^ 1 



[ITU'- 



(r - 2M) 



(r-2M)4 1 )+r^ (0) / t ^ 



(1) 



(14) 



The second order correction to the radial velocity is found 
from the terms at 0(e 2 ) of Eq. (|rU|) : 



[LVJ 



r r 3 w d) ad 



3r 3 o; 2 0) 



(r - 2M) 
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2r 3 (r - 2AfV«4 1 ^ (0) 
lir 2 V {1) Vluj {X) (r - 3M) + Aa/ (1) 



r^ 2 - (r - 2M)1 / t (2) 



(r - 2M)/W 



r 3 w 2 



(0) 



- (r - 2M) | (15) 



As we shall see below, ui^ 2 ' does not contribute at the 
order to which we are solving the EOM. The explicit 
perturbative solution of the EOM is listed in Appendix 
B. 

We remark that we express all quantities here as func- 
tions of r. For other purposes, e.g., for analysis of de- 
tected signals, it is frequently more convenient to express 
quantities as functions of u. It is easy to translate our 
expressions to functions of u by noting that duj — u>' dr. 



involving / t /r in Eq. (|B4|) . The ratio of these two 
terms is 



B t f™ [M/n)B t fW 
7c := — — ~ a t - 



Atrft \fr ' A-trft ^ fr ^ 



{M/n)B t fP 
~ at A f (D • (16) 

Next, for quasi-circular orbits at r > M, ft/fr ~ 
(M/r) 1 / 2 , such that 



7c ~ a t 



{M/n)B t (M 



Ar V r 



a t r-6M (M 
Yr-3M \~r 



1/2 



5/2 



(17) 



B. Approximating the solution and its validity 
regime 

Notice that each term inside the curly brackets in 
Eq. i|15|l is a t 0(/i 3 ). These terms come in two kinds: 
first, there are terms which are quadratic in (or their 
gradients), and second, there are terms which are lin- 

(2) 

ear in fjj, ■ The latter, of course are as yet unknown. 
Their derivation requires second-order perturbation the- 
ory, and the extension of regularization techniques to 
that order. One should therefore view Eq. (|15fl as a 
formal expression. However, we propose that the ex- 
pression for is dominated by the terms quadratic 
in fji, such that neglecting the terms involving /^ 2 ' in- 
troduces only a small error: Consider the t, r components 
of the four-acceleration and the four-force. Expand the 
four-acceleration as a = aW + a^ 2 ' + ... where a^ n ' is 
at order e". Then obviously oc aW , with some pro- 
portionality constant (with dimensions of 1/M) which we 
expect to be neither very large nor very small. (Recall 
that we solve perturbatively, such that ~ jj,/M 2 . In 
this case a*- 2 -* <C a^\) The force = fj,a,u can be 



written as cx fiai, 1 ' . Substituting a\j L> — fv ±! /ft we 
then find that cx /i -1 /^ . As fv has the same 
dimensions as fu (namely, it is dimensionless in a nor- 

malized basis), we expect that /„ = a v M^x fb , (no 
summation over repeated indices implied) where a v is an 
(unknown) dimensionless function of r/M. Next consider 

the ip component. We expect = —f^/uj. Substi- 
tuting Kepler's law for u, and repeating our arguments 

above, we find that — a v (M /r) 3 / 2 yT 1 f$' , where 
again a v is a dimensionless function of r/M. 

Consider now the terms in Eq. (|15|) . The coefficients of 

(2) 

the terms proportional to fp\ are typically much smaller 
than the coefficients of the terms quadratic in fjp. For 

(2) 

example, compare the term involving / t and the term 



P (i) 



Here, B t = 2r(r - 3M)/\p(r - 6M)] and A tr = 4r 3 (r - 
3M) 2 /[fj, 2 M(r - 6M) 2 ] are the coefficients of the terms 

proportional to and f^f^K respectively, in V^ 2 \ 
Demanding that 7c <C 1 introduces then retrictions on 
the function a t . Specifically, for r>M our assumption 
is justified if a t increases as a function of r/M slower 
than (r/M) 5 / 2 . (Even if at violates this condition, our 
assumption can still be valid in a neighborhood of a point 
r at which it is valid. If the condition is satisfied, our as- 
sumptions are valid globally.) Presently, as was already 
discussed above, we have no knowledge of the functions 
a^. We proceed by assuming that a M are such that in- 

(2) 

deed the terms proportional to fp\ are small compared 
with the terms quadratic in in Ea. (|15|l . Similarly, we 
can introduce analogous constraints on a r and a v such 

that any term linear in in Eq. 1)15(1 would indeed be 

much smaller than any term quadratic in ffp . Although 
our discussion above is limited to r > M, we assume the 
applicability of our conclusions for all r > 6M. We there- 
fore neglect all the terms in Eq. (|15f) which are linear in 

f (2) 

We emphasize that even if our assumption that all the 
terms in Eq. (|15|) which are linear in are negligible 

compared with terms which are quadratic in /jp turns 
out to be incorrect, our main conclusion in this paper is 
still likely to be relevant for at least a portion of the orbit 
(see below). 

The solution we find for the EOM is a perturbative 
solution about a circular orbit. So long as the orbit does 
not get too far from circularity we therefore expect our 
solution to be valid. However, when the particle arrives 
at the innermost stable circular orbit (ISCO) at r = 6M 
the orbit changes from an adiabatic, quasi-circular orbit 
to a plunge. (We do not consider here the change in 
the ISCO itself under radiation reaction and spin-orbit 
effects.) When that happens the radial velocity can no 
longer be considered as small, and the perturbative ap- 
proach to the solution (i.e., expanding the solution about 
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a circular orbit) breaks down. To find where our pertur- 
bative approach loses its validity let us compare the mag- 
nitude of V^ 2 ' with the magnitude of V^ 1 ' . We consider 
the perturbative approach as valid only if is smaller 
than V^K Specifically, let us introduce the condition 
that \V^/V^ \ < 7 for the perturbative approach to be 
valid, where 7 is a positive constant smaller than unity. 
(In practice, we arbitrarily fix 7 = 0.2.) For a typical 

term in we take the term proportional to ft fS^ in 
Eq. i|B4(l . and for a typical term in we take the term 

proportional to in Eq. (|fj3() . Their ratio is 



y(i) 



2 r 2 (r-3M) 



6M 



f(i) 



< 



(18) 



As fr 1 ^ is bounded, clearly this inequality is violated at 
some value of r > 6M. However, as fr is at (3(/x 2 ), the 
domain of validity of the perturbative approach extends 
to smaller values of r the smaller /1. 



IV. COMPUTATION OF THE ORBIT AND THE 
WAVEFORM 

In the previous section we found V(r) to order ^t 2 and 
u>(r) to order \i for quasi-circular equatorial orbits. Next, 
we find the orbit by computing 



i(r) 



dr 

W) 



and </?( r ) 



tu(r) dr 
V(f 



(19) 

where the integrands are evaluated using the perturba- 
tive solution for V(r) and ui{r). In practice, we solve 
these integrals using a fourth-order Runge-Kutta inte- 
grator. This yields the triad r, t(r), (p(r), which we invert 
to t, rit), y>(t), which is just the required orbit. 

Next, we compute the number of cycles Af cyc spent in 
a logarithmic interval of frequency /, dAf cyc / din/, by 
noting that dAf cyc / d\nf = u> /(ttuj). Re-expressing the 
last equation as dAf cyc / d\nf = lo 2 / (ttlo'V) and using the 
perturbative solution for V(t') and for u>(r), we find that 



dJ\f cyc 



2 

'3tt 



M 



2^r- + - 



2o/« 



+ 0(/i 2 ) (20) 



Finally, the total change in the number of cycles can 
be found by integrating 



AJVc, 



uj din/ 



■ df . 



(21) 



V. MODELING THE SELF FORCE 

A. Fitting the luminosity in gravitational waves to 
a smooth function 

The local gravitational self force has not been com- 
puted yet for a point particle in motion around a 
Schwarzschild black hole, not even for quasi-circular, 
equatorial orbits. (An exception is a radial plunge, for 
which it has been calculated 01-) However, the ra- 
diation reaction for quasi-circular, equatorial orbits in 
Schwarzschild (and also for other types of orbits, includ- 
ing in Kerr) was calculated using balance arguments, and 
the results for the rate of change of the constants of mo- 
tion can be translated into the local self force. As the 
local self force is not available to us, we shall proceed 
by using results obtained from the non-local approach. 
That approach can give us the dissipative part of the self 
force, but not the conservative part, which is not encoded 
in the gravitational-wave luminosity. We will therefore 
estimate the conservative piece of the self force by anal- 
ogy the the case of scalar field self interaction. Although 
we do not expect our estimate to be accurate, we nev- 
ertheless hope that it can still provide us with a crude 
order-of-magnitude estimate for the actual conservative 
piece of the self force. 

In practice, we take the results for the luminosity of 
gravitational waves (i.e., the flux to infinity) for a test 
mass in circular orbit around a Schwarzschild black hole 
PH . These results are given in the form of a table, which 
gives dE / dt(r). For the purpose of integrating a slowly- 
evolving orbit, we need a smooth function. We can ob- 
tain such a function by fitting the data in the table to 
a smooth function. The most natural functional form 
to use is that of the post-Newtonian (PN) expansion. 
While the PN expansion converges very rapidly at large 
distances, it does so only very poorly close to the ISCO, 
where the motion is very relativistic. On the other hand, 
a general polynomial fit converges very rapidly at small 
distances, but poorly at large distances. We find that the 
two asymptotic expansions match at an overlap region, 
such that a matched asymptotic expansion method can 
be very useful. 

In practice, we use the 3.5PN expansion for a test par- 
ticle at large distances. The flux of energy to infinity is 
given by (l3lll4| 



d£ 3 - 5PN 
dt 



32 n 2 M 3 



1 



1247 



8191 
"672 

856 
105 



-irv 



ln(16?; 2 ) 



336 
6643739519 
69854400 ~ 
16285 



Airv 6 



44711 



9072 



1712 

T05 



■7c + 



16 



504 



-irv 



(22) 



where v = (Muj) 1 ^ 3 is the orbital velocity, and where 7 e 
is Euler's constant. 

At small distances we fit the luminosity to a polyno- 
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FIG. 2: Luminosity in gravitational waves. Upper panel: The 
temporal component of the self force [Eq. I|260 ] as a function 
of r/M. The numerical data are represented by CPs, and the 
fitted curve by the solid line. Lower panel: The relative error 
in the determination of the luminosity F using the 3.5PN ex- 
pansion [□, F := ( dE/ d£) 3 ' 5PN ] and the power-law expansion 
[o, F := (dE/dt) PL ]. 



mial of the form 



dE 
~~dt 



PL 



32 M 2 M3 » 



(23) 



fc=0 



with coefficients Ck given in Appendix IU1 

Figure [5] shows the temporal component of the self 
force and the relative errors in its determination in terms 
of both methods. As expected, the relative error in the 
3.5PN expansion is very small at large distances, but 
grows rapidly at small distances. Similarly, the power- 
law approximation yields highly accurate results at small 
distances, but is inaccurate at large distances. There 
is an overlap region, for 12 < r/M < 24, where the 
two expansions have comparable accuracy. Specifically, 
the 3.5PN expansion yields at r = 12M a relative er- 
ror of 1 - (dE/dt) 3 - 5PN /(dE/dt) « 4 x 1CT 4 , and the 
power-law expansion yields at r = 24M a relative error 
of 1 - (dE/ dt) PL /{dE/ dt) w 1.1 x 1(T 3 . Allowing our 
accuracy in the determination of the luminosity to be 1 
part in 10 3 , we evaluate the self force at large distances 
(r/M > 24) using the 3.5PN expansion, and at small dis- 
tances (r/M < 12) using the power-law expansion. For 
the region 12 ^ r/M ^ 24 we use a sliding average of the 
two methods, i.e., we simply take 



dE' 
dt 



dE PL 24:M-r dE 



3.5PN, 



12M 



dt 



12M 



dt 



12M 



(24) 



This guarantees that our fit to the luminosity, and to the 
temporal component of the self force, will nowhere be 
worse than our tolerance. 



B. Error analysis 

As we already discussed in the preceding subsection, 
the error in the determination of the luminosity in gravi- 
tational waves for a test mass is no worse that 10 -3 . The 
other sources of errors which affect our determination of 
the ^-component of the self force are finite mass effects 
(i.e., the effect of a finite (t/M), and the flux of energy 
through the black hole's event horizon. 

The total rate of change of energy of the particle is the 
sum of the luminosity in gravitational waves (i.e, the flux 
of energy to infinity) and the flux in energy through the 
event horizon of the black hole. Specifically, 



dE TOT _ dE°° 
dt ~ dt 



dE EK 
dt 



(25) 



As ( dE/ dt) Eii m av 8 (dE/ dt)°°, with a being a constant 
of order unity 15], the flux of energy though the event 
horizon where it is the greatest (when the particle is at 
r — 6M) is smaller than the flux to infinity by a factor 
of 8 x 10~ 4 . Consequently, neglecting the absorption by 
the black hole introduces an error in the determination 
of the ^-component of the self force which is compatible 
with our tolerance. 

The other source for errors in the determination of 
the ^-component of the self force are finite-mass effects. 
The finite-mass effects are easy to estimate for large dis- 
tances, where they are fully known at the 3.5PN level. 
(Recall however the ambiguity in determination of the 
finite-mass effect of the 3PN terms.) The error intro- 
duced by neglecting the finite-mass effects in the 3.5PN 
expansion is (a few) x fi/M. The error will be compat- 
ible with our tolerance if we take the mass ratio to be 
\i/M < (a few) x 10 -4 . We shall therefore restrict the 
systems we study here to compatible mass ratios. We 
cannot estimate finite-mass effects close to the ISCO, as 
the luminosities available to us were obtained from a lin- 
earized analysis. We are encouraged, however, that at 
all PN orders (up to 3.5PN) the errors introduced by 
neglecting finite-mass effects are comparable. If that be- 
havior persists at all PN orders, we are guaranteed that 
our approximation is valid also near the ISCO. Although 
we do not know whether this is indeed the case, we hope 
that our approximation remains qualitatively valid also 
near the ISCO, and that it does not introduce errors sig- 
nificantly bigger than our tolerance. 



Determination of the self force 



by 



The ^-component of the self force is easily determined 



SF 



dt dE 
dr dt 



(26) 



In Eq. (|26|l dE/ dt is the flux of energy to infinity. 
The rate of change of the particle's energy then is 
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— dE/dt. Notice, that we only need to determine / t , 
such that in Eq. (|26l) the factor dtj dr can be evaluated 
for geodesic motion, i.e., dtj dr = 1/(1 — SM/r) 1 / 2 . The 
(^-component of the self force can be determined from its 
proportionality to the i-component. 

The r-component of the self force cannot be found by 
using non-local methods. Instead, it needs to be eval- 
uated using a fully local calculation. As the results of 
such a calculation are as yet unavailable to us, we shall 
instead use a crude order-of-magnitude estimate for f~ F . 
Because of the smallness of the effects of this compo- 
nent, and because it is a posteriori found to be much 
smaller than the magnitude of the r-component of the 
spin force, we do not expect our crude estimate to be 
problematic. Note, that once is calculated, it can 
easily be used to replace our estimate here. Specifically, 
we estimate the radial component of the self force by us- 
ing the known results for the radial self force on a scalar 
char ge q in circular orbit around a Schwarzschild black 
hole |l6lll7| . In order to do that estimate, recall that for 

r > M, grav/f/scalar/f « fWg) 2 (M/r) 0, G3 . We 

next assume that an analogous scaling is satisfied also by 

(2) 

ffi, (which is entirely conjectural, although not implau- 
sible). Consequently, because /f F u" = 0, it is reason- 
able that the same proportionality is satisfied also by the 
r-components, fl X \ For strict circular orbits the radial 
velocity vanishes, such that no restrictions on can be 
applied, but for quasi-circular orbits the orthogonality of 
the force and the orbit imply that it is not implausible 
that such a proportionality is satisfied. In view of this 
argument, we assume 

fSF _ v fSF / 97A 

gravj t ~ " 9 A scalar,/ r \" 1 ) 

qz r 

Although we expect f3 = 96/5, in practice we parametrize 
our results with (3. Although our arguments may apply 
to the far field region only, we extrapolate the scaling re- 
lation l|27(l for all values of r. This choice, in the absence 
of numerical data, appears to us to be reasonable, at least 
as an approximation for the actual radial component of 
the self force. 



VI. CALCULATED ORBITS AND WAVEFORMS 

In order to compare between different evolutions we 
use dhf cyc / dhxf, the number of cycles Af cyc that the 
system spends in a logarithmic interval of frequency /, 
as given in Eq. I|20|) . The total number of cycles that the 
system undergoes between two values of r, AAf cyc (this, 
of course, can be translated to the number of cycles be- 
tween two values of the frequency /), is found by using 
Eq. H21f> . First, we discuss the contribution to AJ\f cyc 
from the spin-orbit coupling. Figure |21 shows the differ- 
ence between AAf cyc which was obtained for s = 1 and a 
number of values of s, for a system with \x = 5 X 10 M 
that starts decaying at r = 10M. The maximum effect 




6 6.5 7 7.5 8 8.5 9 9.5 10 



r/M 

FIG. 3: The difference between AA/" C yc which was obtained 
for s = 1 and a number of values of s, as a function of r/M. 
The data are for a system with /j = 5x 10~ 4 M that starts 
decaying at r = 10M. Shown are the results for s = — 1 (solid 
curve), s = —0.5 (dashed curve), s — (dotted curve), and 
s = 0.5 (dashed-dotted curve). 



is obtained when s = — 1, where it is just above 1 cycle 
at r = 6M. This effect is at 0(/x°), i.e., it is independent 
of the mass ratio p|. That is, ignoring the spin force 
would result in a maximal error of about a full cycle, 
which could reduce the correlation integral of the signal 
with a theoretical template considerably. (Notice that 
the total number of cycles that the system undergoes is 
A/" C yc ~ 2.6 x 10 3 , and that J\f cyc oc fi^ 1 .) 

The waveforms, which we present using the "restricted 
waveform" approximation [lg, are displays in Fig. IVII 
(Our calculation of the AA/" cyc is independent of the as- 
sumption of the "restricted waveform" approximation.) 
Figure IVII indeed shows that the maximal ambiguity in 
phase due to ignoring the spin-orbit coupling is about 
two wavelengths. 

In order to examine how important that self force may 
be we fix s = 1, and compute 5[AJ\f cyc ], the difference 
in the total number of cycles for the case in which we 
include the self force in the calculations of the orbital 
evolution, and the case in which we do not. The results, 
for j3 = 1, are presented in Fig. [5] We find that the 
difference in the total number of cycles is 2 x 10 -3 /? at 
r = 6M. Even for j3 ss 20 the difference in the total 
number of cycles is only 4 x 10 -2 , which is very small 
indeed. This effect is at O(fjP) @. Figure displays a 
fraction of the last oscillation in the waveforms with and 
without the inclusion of the self force (for s = 1) before 
r = 6M. 

Most importantly, we can make the waveform without 
the inclusion of the self force coincide almost exactly with 
the waveform with its inclusion, if we modify the value of 
s by a small number which is comparable to the difference 
in the number of cycles for the cases with and without 
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1.816 1.8165 1.817 1.8175 1.818 1.8185 

t/M X10 5 

FIG. 4: The waveform - h as a function of t/M. The data are 
shown for /i = 5 x 10 _4 M. We show the last few oscillations 
before r = 6M. Plotted are the waveforms for two values of 
s: s = 1 (solid curve) and s = — 1 (dashed curve). The two 
curves were in phase at t = (r = 10M). 




6 6.5 7 7.5 8 8.5 9 9.5 10 



r/M 

FIG. 5: 5[A7V C yc] — the difference between AJ\f cyc with the 
inclusion of the self force, and ATVcyc without its inclusion — 
as a function of r/M. The data shown are for n = 5 x 10" 4 M, 
s = 1, and the system starts decaying at r — 10M. 



the inclusion of the self force. This situation is shown 
in Fig. [7] which displays two waveforms: one waveform 
is the same as the waveform without including the self 
force effect in Fig. [5] The other waveform is a waveform 
with the self force effect, but with a corrected value of 
the spin s. The two waveforms overlap almost exactly 
over the entire wave train, and in particular during the 
last fraction of the orbit before r = QM . 

Therefore, we propose that neglecting the self force will 
only introduce a small error in the determination of the 
spin rate of the companion, at the order of <5[AA/^ yc ]. 



0.17 - 




0.07 - 

-8 -7 -6 -5 -4 -3 -2 -1 

t/M 

FIG. 6: The waveforms for fj, = 5 x 10~ 4 M and s = 1. Solid 
curve: including the effects of the self force. Dashed curve: 
neglecting the effect of the self force. Here, we used j3 = 1, 
and the two waveforms were in phase at r — 10. Shown is a 
fraction of the last oscillation before r — 6M. Wr set t = at 
r = QM. 



0.17 - 




0.07 - 

-8 -7 -6 -5 -4 -3 -2 -1 

t/M 

FIG. 7: The waveforms for y. = 5 x I0 _4 M. Solid curve: 
excluding the effects of the self force for s = 1. Dashed curve: 
including the effect of the self force for s = 0.997. Here, we 
used f3 — 1, and the two waveforms were in phase at r — 10. 
Shown is a fraction of the last oscillation before r = 6M. The 
two waveforms are indistinguishable on this scale. 
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APPENDIX A: EXPRESSIONS FOR /f F AND ff pln 



The Papapetrou equations are then given by 



In this Appendix we write explicitly the expressions for 
/gp and for / s ^, in which appear in Eq. Specifically, 

/ S " F = fi 2 k a ^ 5 lim f V 5 G wv ^(r);^(r')] 



x u? \t')v?' '(r')dr' 



(Al) 



where 



2 y y 2 



+ -g^u a u 5 + \g aB a 01 



(A2) 



and G / 3 7( 3' 7 '[z' 1 (t}; ^(t')] is the two-point retarded 
Green's function |g. 

The spin force is given (at the pole-dipole approx- 
imation, i.e., taking into consideration only the mass 
monopole and the spin dipole, neglecting higher mul- 
tipoles, such as the tidal coupling, which is a mass 
quadrupole effect) by the Papapetrou equations 



and 



Dp a 1 
dr ~ 2 b H 



DS al3 n r„ DS^p 



dr 



dr 



(A3) 



(A4) 



Du a 
dr 



Xfipa jDaiy 



and 



DS a _ l)W /wc2i 



dr 



dr 



= + O(S z ) 



(A8) 



(A9) 



which means that to leading order in the spin, the Pauli- 
Lubanski spin covector is parallel transported along u a . 
(The spin covector is always Fermi- Walker transported.) 
In these equations we neglect all terms at order S 2 or 
higher. From Eq. I|A8(1 we identify the spin force as 



1 



f a — _A^P< J R al/ „ „ c 
./spin — 2 \ tl a v a <T->p 



(A10) 



APPENDIX B: EXPLICIT EXPRESSIONS FOR 
SCHWARZSCHILD 

We list in this Appendix the explicit expressions for 
the solution of the perturbative equations of motion. We 
remark that the following explicit solution appears dif- 
ferent from the solutions of Ref. [j| because of the differ- 
ent perturbative expansion of the angular velocity that 
is used there. 



where p^ — (iu a — up d ^ t is the particle's four- 
momentum, S a/3 is the skew-symmetric spin tensor of 
the particle, which is given by 



3Af 



2 J u(Mr) 1 / 2 



(i) 



(Bl) 



S aP = 2 / (x [a - z^ a )T^ dE. 



7 ' 



(A5) 



where £ is an arbitrary spacelike hypersurface. Notice 
that the RHS of Eq. is just 2p^ a u^. 

As the mass of the particle /i <C M, and as the spin of 
the particle is consequently also small, we shall approx- 
imate the Papapetrou equations by considering the spin 
force to leading order, that is to order /i 2 We next 
introduce the Mathisson-Pirani spin supplementary con- 
dition (which physically identifies the particle's center of 
mass) S a @up = and the Pauli-Lubanski spin covector 
S a , defined by 



S a 



1 



(A6) 



(In fact, to leading order in the spin the Mathisson- 
Pirani spin supplementary condition is identical with the 
Tulczyjew-Dixon condition S a @pf3 = 0. That is, the dif- 
ference between the two is cubic in the spin of the par- 
ticle, and hence negligible.) Notice, that the inverse of 
Eq. (2H> is given by 



5 



a/3 



U-ySg 



(A7) 



V(Afr) 1 /2 



r + 3M (1 



fP + 2(r - 3M )/«' 



(B2) 



2r r-3M 
JiMr- 6M 



r J 



fr+Mfi 



(i) 



(B3) 



(2) 



r(r - 3M) 
fi 2 M 2 (r-6M) 2 



X (r - 3M) + l—j ' 4 1} /i 1} (5r - 6M)(r - 2M) 

X (r - 3M) + 2M fP r 2 {r - 2M)(r - 3M) 
+ 4M/ t (1) /i 1 V 2 (r - Ml) + 2^M 2 fP(r - 6M) 
' M 



ff\r-2M){r-m) 



(B4) 
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APPENDIX C: THE EXPANSION 
COEFFICIENTS C k 

We fit the luminosity in gravitational waves (data 
taken from Ref. to a function of the form l|23|) with 
n = 6. The number of terms we sum over was determined 
by the accuracy of the approximation. The coefficients 
were determined to equal (not all the figures shown are 
significant) 



Ci 
C 2 

c 3 
c 4 
c 5 

Cc, 



= -0.9253608775 
= 0.1300173368 
= -9.256723735 x 10~ 
= 3.455531402 x 10~ 4 
= -6.363613595 x 10" 
= 4.521868774 x 10~ 8 



(CI) 



Co = 3.604707682 
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